Weak Vàlues, Quantum Trajectòries, and the Stony-Brook Cavity QED experiment 
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Weak vàlues as introduced by Aharonov, Albert and Vaidman (AAV) are ensemble average vàlues 
for the results of weak measurements. They are interesting when the ensemble is preselected on a 
particular initial state and postselected on a particular final measurement result. I show that weak 
vàlues arise naturally in quantum òptics, as weak measurements occur whenever an open system is 
monitored (as by a photodetector) . I use quantum trajectory theory to derive a generalization of 
AAV's formula to include (a) mixed initial conditions, (b) nonunitary evolution, (c) a generalized 
(non-projective) final measurement, and (d) a non-back-action-evading weak measurement. I apply 
this theory to the recent Stony-Brook cavity QED experiment demonstrating wave-particle duality 
[G.T. Foster, L.A. Orozco, H.M. Castro-Beltran, and H.J. Carmichael, Phys. Rev. Lett. 85, 3149 
(2000)]. I show that the "fractional" correlation function measured in that experiment can be recast 
as a weak value in a form as simple as that introduced by AAV. 
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I. INTRODUCTION 



The concept of weak vàlues in quantum mechanics was 
formulated by Aharonov, Albert and Vaidman (AAV) [ 1 1 , 
using the earlier two-wave-function formalism of Ref . [ 2 1 . 
A weak value is the ensemble average value of a weak 
measurement result. A weak measurement is one that 
minimally disrupts the system, while consequently yield- 
ing a minimal amount of information about the observ- 
able measured. For a given initial system state, the 
ensemble average of weak measurement results is the 
same as for strong (i.e. projective) measurement results. 
Where weak measurements are interesting is when a fi- 
nal as well as an initial state is specified. Here the final 
state is the result of a second measurement (a strong 
one), so that the ensemble average is taken over a post- 
selected ensemble, in which the desired result for the final 
measurement was obtained. 

Let the desired initial system state be and let i/j 
denote the event that this is successfully prepared. Let 
the desired final system state be \4>), and let <f> denote 
the event that the appropriate result is obtained by the 
final projective measurement. Let the observable to be 
measured at an intermediate time have operator X, and 
denote the result of the measurement X n . Here the na- 
ture n of the measurement could be strong (n = s) or 
weak (n = w). In the former case, which corresponds to 
a projective measurement, X s reproduces the statistics 
of X. In the latter case, X w does not, and indeed is not 
confined to the spectrum of X. We do require however 
that in all cases X n be an unbiased estimator, in the 
sense that 



(X n ), = E[X n \i;} = 



(1.1) 



Here E denotes expectation value, and the vertical line 
denotes "given" . The postselected value we desire is 



,{X n )^ = E[X n \^,<p] = 5>Pr[A„ = x\cf>^] 

X 

= E x xPr[c^\X n = x, jj]Pv[X n = a#] 
V Pr[0,A„ = y|^] 



(1.2) 
(1.3) 



Here Pr[ü] denotes the probability for event H . 

For simplicity, assume for the moment that there is no 
evolution between measurements. Then for a strong in- 
termediate measurement, the postselected average value 
of X s is 

/y v T, x ^[Ms=x]Pr[X s =x\^] 



£ y Pr[0|A s =y]Pr[A s = y|^] 

EJ(*)MW)I 2 



EJ(*)I 2 IW)I 2 



(1.5) 



Here \x) denotes an eigenstate of X. In the final result 
there is no direct connection between the initial state and 
the final state. The intermediate strong measurement has 
destroyed that connection, as it enables the dr opp ing of 
the conditional \ijj from Pr[0|A s = x,ip] in Eq. (1.3). Ba- 
sically there are no surprises here. 

Contrast this result with that from a weak measure- 
ment, as derived by AAV: 



Re 



(1.6) 



Aharonov and Vaidman derive this in detail in Ref. ||, 
and it is also derived here, in Sec. II, from a somewhat dif- 
ferent perspective. This formula appears to have a muc h 
more of a mysterious "quantum" nature than Eq. (|l.5|), 
as it involves the initial and final states linearly (rather 
than bilinearly) on the numerator and denominator. As a 
consequence, weak vàlues can have bizarre properties, es- 
pecially when the outeome <f> used for postselection is im- 
probable ||. For example, the weak value ^{Xw)^ may 
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lie outside the range of eigenvalues for X . Counter- 
factual paradoxes kjj can also be expressed and resolved 
in this language ||Q . 

It might be thought that weak vàlues would be mainly 
of theoretical interest, because arbitrarily weak measure- 
ments would not arise naturally in experiments 0. But 
in fact arbitrarily weak measurements arise all the time, 
whenever a system is monitored (i.e. measured con- 
tinuously in time). The measurement theory for mon- 
itored systems was developed by mathematical physi- 
cists [^-[ÏÏ) • It was subsequently rediscovered by workers 
in quantum òptics ]Ï2]-|Ï4| when it was found necessary 
to describe the conditioning of systems on an observed 
photocurrent. In quantum òptics, the theory of continu- 
ous monitoring, especially for homodyne and hctcrodyne 
measurement, has come to be known as quantum trajec- 
tory theory [[14 15 1. 

In this paper I show in Sec. II how the AAV formula 
can be derived using the approach of quantum trajectory 
theory. I derive a simple generalization arising from re- 
moving the restriction that the weak measurement of X 
be a quantum non-demolition measurement. In Sec. III 
I derive a more complete generalization which allows for 
(a) mixed initial conditions; (b) non-unitary evolution 
between measurements; and (c) a non-projective final 
measurement. Within this context, the AAV formula can 
be seen to be a special case of the correlation functions 
that have been used in quantum òptics for a long time 
p6| . Finally, in Sec. IV I show that the correlation func- 
tion measured in the recent cavity-QED experiment at 
Stony-Brook |Ï7j can be reformulated as a weak value 
almost exactly as originally proposed by AAV. 



II. WEAK VÀLUES: THE SPECIAL THEORY 

Monitoring a quantum system can be treated by con- 
sidering a scqueiice of discretc weak measurements, each 
taking time St, and letting St — > dt. To obtain a sensi- 
ble limit, the strength of each measurement must scale 
as St. By this I mean that the post-measurement sys- 
tem state should, on average, be different from the pre- 
measurement state by an amount of order 5t. 

To describe weak measurements, we require gener- 
alized (non-projective) quantum measurement theory 
pi| pof . Let the measurement result be the random vari- 
able X n as above. In the case of efficient measurements, 
the measurement can be described in terms of a set of 
measurement operators {M x } x , each associated with a 
result X n — x. The post-measurement conditioned state 
is given by 



By conservation of probability, the measurement opera- 
tors must obey a completeness relation 



Px{t + St) 



M x p{t)Ml 



Pr[X„ = x] ' 
where the probability for the result is 

Pr[X n =x]=Tr[M x p(t)M x ]. 



(2.1) 



(2.2) 



£4 = ï, 



(2.3) 



where É x — M\M X is known as the •probability operator 
or effect for the result x. 

As above, we require the measurement of X to be un- 
biased, so that (X n (t)) p ^ = Tr[p(t)X]. This is satisfied 
by the following measurement operator, which also has a 
suitably weak effect on the system (^ÏJ 

M x = (2ir/5ty 1,i ex-p(-x 2 8t/2) [l + St(xc - c f c/2)] . 

(2.4) 

Here c is an operator such that 

c + c)=X. (2.5) 

For convenience we are measuring time in units such that 
&c St is dimensionless. Note that by allowing c to be non- 
Hermitian we are allowing for non-back-action-evading 
measurements of X, even ignoring any evolution between 
measurements . 

The measurement operators sa tisfy a continuous ver- 
sion of the completeness relation ( |2.3| ) to order (St) 2 : 

J dxM\M x = 1 + 0(St 2 ), (2.6) 

In the limit St — > dt, the correction vanishes and it is not 
difficult to show J2ÏJ that the measurement result X w has 
statistics given by 



X w (t) = Tr[Xp{t)]+Ç{t). 



(2.7) 



Here is a Gaussian white noise term, defined by £(t) 
dW(t)/dt where 



E[dW(t)} = 0, 
dW{t) 2 = dt. 



(2.8) 
(2.9) 



Moreover, the stochastic post-measurement conditioned 
state is given by 

pxjt + dt) = {l + V[c]dt + H[c]dW(t)}p{t). (2.10) 

Here for arbitrary operators A and B 

v\A]È = AèÀ^-\{À^À,È}, (2.ii) 

H[À]È = (À-Tt[ÀB])B + U.c. (2.12) 

Averaging over the measurement results would give 

dp = dW[c]p, (2.13) 

which does not necessarily preserve (V) p ( t ) = Tr[A /?(£)] 
because c need not commute with X . 
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For continuous measurements we wish to consider a 
sequence of measurements as described above. However, 
for the purposes of deriving weak vàlues, we wish to as- 
sume we have just a single such weak measurement, at 
time t. (This unrealistic assumption will be removed in 
Sec. III.) Say the system was prepared at time in state 
IV'(O)) and evolved unitarily to state \tp) at time t. Let 
us use p(t) for \il))(ip\. The infmitely weak (i.e. infinitesi- 
mally strong) measurement is then performed [ p2| . Then 
at time T a final projective measurement is made, and we 
wish to keep only those vàlues of X obtained for which 
the final result corresponds to the state |çí>(T')) . Let \<f>) 
be that state evolved unitarily backwards in time from T 
to just after t. Dropping the time arguments, we have in 
the notation of Sec. I, 

=E[X w \4>,Tp}=J2x?*[Xw=x\<t>,Í>} (2-14) 

X 

= 5^a;Pr[X„=x,^]/Pr[^]. (2.15) 

X 

Here we have used a sum rather than an integral merely 
for convenience. 

Now because the measurement is infmitely weak, it 
does not affect the denominator D = Pr[0|ï/>], which 
equals |(0|í/')| 2 . The numerator meanwhile evaluates to 

N = Y^ x¥ï[(j)\X w = x, i>]Pr[X w = a#] (2.16) 

X 

= Y d ^<t>\px{t + dt)\<j>)Vx[X w =x\^\ (2.17) 

X 

= E[X w (cl>\px w (t + dt)\<P)]. (2.18) 



Now from the above formulae (p.7|)-(2.l·0) it is easy to 
see that this evaluates to 

(0| {p(t)(X) p{t) + H[c]p(t)} \<t>) + O(dt). (2.19) 

Expanding the superoperator TL and ignoring infinitesi- 
mals gives 



(<P\{cp(t)+p(t)c^} 



(2.20) 



Substituting this into Eq. (2.15), and using p{t) — \ïp)(ip\, 
gives 



2Re 



■\c\^)(ip\4>) +H.c. 

\(M)\ 2 



(2.21) 



(2.22) 



In the case of QND measure men ts where c — X /2, this 
reduces to the AAV formula (1.6). Putting in the unitary 
time evolution explicitly, and setting T = t + r, 



,{X w {t)), = 2Rc 



^(T)\Ü(r)cÜ(t)\^(0)) 
{0(T)|{7(T)|V(O)) 



(2.23) 



III. WEAK VÀLUES: THE GENERAL THEORY 

We now generalize the above theory by removing many 
of the assumptions. First, we allow the initial system 
state p(0) to be mixed. Second, we allow arbitrary 
Markovian evolution [pp[ of p between measurements: 



P 



Cp=-i[H,p] + ^2v[c„]. 



(3.1) 



Here H is an Hermitian operator while {c M }^ is a set 
of arbitrary operators (which strictly should be bounded 
Ipijl). It would be natural for one of them to be the op- 
erator c introduced in Sec. II. Then at any time between 
preparation and final measurement, the conditional evo- 
lution of the system would be 



dpx m ( s ) — dsCp + dW(s)7í[c]p. 



(3.2) 



Finally, we allow for the final measurement at time T 
to be described by a positive operator E(T) which is 
not nccessarily a projector \(j>(T)}((j)(T)\. As long as 
1 — E(T) > 0, E(T) is an effect for some final measure- 
ment. 

Using E and p to denote the events of successful prepa- 
ration and successful final measurement, the weak value 
at time t we wish to calculate now is 

E (X w ) p = E[X w \E,p] =J2^4^w = x\E,p} (3.3) 

X 

= xPx[X w - x, E\p}/Pr[E\p\. (3.4) 

X 

As before, the denominator is unaffected by the weak 
measurement and is given by 



D = Tr[E(T)e L·1 p(0)]. 



(3.5) 



The numerator can be manipulated similarly to above to 
yield 



N = E X w (t)Tr[È(T)e CT p Xw (t + dt) 
= Tï[È{T)e CT {ce £t p(0) + H.c.}]. 

Equation ( [2.23 ) thus generalizes to 

Tr[È(T)e CT (ce £ V(0)}] 



: (X w (t)) , = 2Re- 



Tr[E(T)e CT p{0)} 



(3.6) 
(3.7) 



(3.8) 



The above result looks more like a normalized correla- 
tion function, as is familia r in quantum òptics, than the 
unexpected AAV formula (L6|). Indeed, it can be written 
in the Heisenberg picture [23j as 



E(T)c(t) + cï(t)E(T) 



(3.9) 
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with initial conditions supplied by p(0 ) and the bath state 
implied by the ma ster equation (|3.l|). However, we can 
formulate Eq. ( |3.8 |) more like Eq. (1.6) as follows. First 



we denote, in analogy to = U(t)\ip(0)), 
P = e Ct P (0). 



(3.10) 



This is th e u sual forward time evolution by the master 
equation (3.1). Next we define a retrodictive |^5j positive 
operator E according to 



Tr[EA] = Tr[E(T)e L·T A], 



(3.11) 



where A is an arbitrary operator. 

It is easy to verify that this E is E(T) evolved forward 
by time r according to the retrodictive master equation 



dE 
~dt 



-i[H,È] 



E [ è t Éè » 



^,E} 



(3.12) 



Note that this differs slightly from the retrodictive master 
equation in Rei. (^6j. This is because we do not require 
the norm of E to be preserved, as it appears in bot h the 
numerator and denominator of Eq. (3.8). That Eq. ( 3.12j ) 



does generate a positive operator is obvious from the so- 
lution as a "Dyson-like" sum constructed in the manner 
of Eq. (7.24) of Ref. |ÏJ]. Using these definitions, we can 
write the postselected weak value as 



: (X w ) o = 2Re 



Tr[Ecp] 



(3.13) 



Thus retrodiction, quantum trajectòries and weak vàlues 
are all united in this expression. 



IV. WEAK VÀLUES IN THE STONY-BROOK 
CAVITY QED EXPERIMENT 

The preceding section showed that postselected weak 
vàlues can be thought of as a correlation function (3.9). 
In quantum òptics weak measurements of the type de- 
scribed in Sec. II are realized through homodyne detec- 
tion [[Ï4 15 1, where c is proportional to the lowering op- 

The typical correlation 



crator of the radiating system. 
functions thus found are of the form (in the Heisenberg 
picture) 

(: X(t + T)X{t) = (c(t + t)c(í) + c\t + r)c(í) 

+ c){t)c(t + t) + c f (í)c + (f + r)) . 

(4.1) 

Usually they are evaluated at steady state (ss) where 

p(0) = p ss oc nullspace(£). (4.2) 



Note that this is not of the form of Eq. ( |3.9| ), so that 
this correlation function cannot be thought of as a weak 
value in the AAV sense. 



A notable exception to this typicality is the recent cav- 
ity QED experiment conducted at Stony-Brook JÏtJ . This 
involved a damped cavity, very weakly driven, through 
which a beam of resonant àtoms passed. The master 
equation for the system can be approximated as 

p = E[à* — à, p] + 2nD[a\p 

N 

+ E - & fa P] + 2 7±2>fo]p} • (4-3) 



Here àj — \g)j{e\ is the lowering operator for the jth 
atom, and à is the annihilation operator for the field. 
The output from the cavity is split into two beams by an 
85-15 beam splitter, One beam is detected by homodyne 
detection (with nett efficicncy r\h < 85%), and the other 
is detected using a photon counter (with nett efficiency 
7 h < 15%). 

Choosing the slightly odd convention of measuring 
time in units of l/(2Kr/h), the homodyne measurement 
gives a current proportional to 



X w (t) = Tr[p(t)X]+Ç(t), 



(4.4) 



where X = à + c 
adding the term 



It conditions the system state p by 



dW(t)H[à]p 



(4.5) 



to dp. Thus the homodyne measurement gives a con- 
tinuous weak measurement as in Sec. II, with c = à. 
The photon counting also gives a continuous measure- 
ment that is weak in the sense that the average change 
in the conditioned system in time St is of order 5t. How- 
ever, unlike homodyne detection, sometimes the change 
is great. This can be seen from the set of effects {E , Ei} 
describing this measurement: 



Èq = 1 — StKr/cà'à, 
Ei = 5tKï] c à^à 



(4.6) 
(4.7) 



These give the probabilities for detecting zero or one pho- 
ton over a short time St. The latter result reveals a lot 
about the system, and its effect (and measurement oper- 
ator) are not close to unity. 

In the Stony-Brook experiment, the correlation be- 
tween the homodyne current X w (t) and the photon count 
increment at time T = t + r were measured in order to 
demonstrate an aspect of wave-particle duality. In the 
experiment r ranged from negative to positive, and the 
measured correlation was proportional to (in the Heisen- 
berg picture) 



h(t - T) 



(à^a)(T)X(t) 



(4.8) 
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This is a "fractional" (3/2) order correlati on f unction [in 
contrast to the first order function in Eq. (4.1)1. Here we 
are interested only in the case r > in which case the 
correlation function can be written as [E7| 



H-r) 



Tr[£^i(T)e £T (àp ss + p ss à 
Tr[É!(T)p ss ] 



(4.9) 



Defining È\ to be Éi(T) retrodicted (as in Sec. III) from 
time T to time t, we can rewrite this correlation function 
as 



Bi (X ro > poB =2Re 



TrfE'ipss 



(4.10) 



It is now clear that the correlation function measured 
in Ref. (ÏJ] is a weak value, preselected by the system 
being in its stationary state, and postselected on the the 
final measurement result E\ (a photon detected at time 
T) . Actually, since what was measured was a correlation 
function (i.e. a function of r), this experiment showed 
the dynamics of a weak value over time. In this case, the 
spectrum of X is the real line, so there is no chance of ob- 
serving a weak value outside of this range. Nevertheless, 
the weak vàlues measured in the experiment were, for a 
range of times r, very far away from the stationary aver- 
age value of X. This lead to the violation of various clas- 
sical inequalities |^,[Ï7). In hindsight, the strangeness of 
the weak vàlues in this experiment is not surprising, since 
the condition cited by Aharonov and Vaidman || is ful- 
filled. That is, the postselection is done on a rare event 
(the detection of one photon rather than zero photons). 

There are a number of features of the above s ystem 
that let us transform the generalized weak value ( 4.10| ) 
into a form almost identical to the special form ( 2.22| ) 
derived in Sec. II. Because the driving is very weak 
(ejC k,j±_), the "jump" terms in the master equation 



(4.3) have a very small effect. This means that it is pos- 
sible to approximate the Liouvillian as 



C = H[-iH tB ] ) 



(4.11) 



where 



iH c ff = e(à^ — à) — c 



N 



f à) - 7x<7 J ta J · 



(4.12) 



Ignoring normalization, this generates nonunitary evolu- 
tion according to the operator 



N(t) = exp(-iH eS t). 



(4.13) 



The stationary solution is, to a very good approximation, 
pure HI 



Pss = |ï/>ss)(-0ss|, 



(4.14) 



where |-0ss) is the eigenstate of —iH g with eigenvalue 
having the largest real part. 

These results mean that we can approximate Eq. ( 4. 10| ) 

by 



Ei (X w )^=2Rc 



ss / 



(4.15) 



Again because of the weak excitation, \i/j bs ) is, to a good 
approximation, given by 



|^ ss ) = |0) + a|l)+/3|l'}. 



(4.16) 



Here a and (3 scale with e. The state |0) is the state with 
no excitations, |1) that with one photon but all àtoms 
in the ground state, and |1') a state with one collec- 
tive atòmic excitation but no photons. The same sort of 
expansion holds for à\ip ss )/ \J {à^à) ss (although the pres- 
ence of a nonzero |1) term here indicates that \ipss) does 
have higher order terms). Furthermore, to leading or- 
der, N(t) preserves this expansion (apart from normal- 
ization). For these reasons, we can replace à'à, which is 
really à^à <g) l a toms, by |1)(1|. That is, postselection on 
the detection of a photon is equivalent to postselection 
on projection into the state |1). The final result thus is 



?1 (x w ) 



2Re 



(i|iv(T)a|v> BB ) 

(1\N(t)\^ s ) 



(4.17) 



This is completely analogous Eq. (F6) when we recognize 
ÍV^(t)|1) as the unnormalized retrodicted final state. 



When \tp ss ) and N(t) are substituted into Eq. (4.17), 
and only the terms of leading order in the driving e are 
kept, the result is 



3i ( X w} % 



{ifj as \X\ip ss ) 



1 + Ce 



-IjT 



cos ílr 



ü 



sin Qt 



(4.18) 



where rj = (k + jj_)/2 and where Ç and Í2 are functions of 
k, 7xj and the coupling coefficients {gj}j- This is exactly 
the same as the analytical result in Ref. ]Ï7fl , where this 
result wa s derived from considering the correlation func- 
tion (i.è) with T < t. That is, with the photodetection 
preceding the homodyne measurement. In that case, the 
correlation function is 



h(-r)=h(\r\) 



Ti[Xe c ^àp ss à^} 
Tr[à/o ss at] 



(4.19) 



In Ref. [Tt]] the symmetry of h(r) was established by 
assuming Gaussian fluctuations. The analysis in this pa- 
per shows that the symmetry (to leading order) of h(r) 
can be demonstrated explicitly by calculating both cases. 
The case of negative r is found simply from the usual pre- 
dictive quantum mechanics with wavefunction collapse 
p9[ , and the case of positive r equally simply by retrod- 
ictive quantum mechanics with weak vàlues, as here. 

Finally, it is important to point out that the above 
form for the weak value as a function of time (4. 18]) 
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matches well with the experiment [[L7|. This is so even 
though the experimental conditions are far from ideal, 
with fluctuations in the number of àtoms, transverse mo- 
tion of the àtoms, imperfect beam alignment and so on. 
This shows that a weak value is not as fragile as one might 
have thought, and does give quantitative predictions for 
the behaviour of the system. 



V. DISCUSSION 

In this paper I have combined quantum trajectory the- 
ory, and the concepts of weak vàlues and retrodiction in 
order to show the relation between weak vàlues and cer- 
tain correlation functions as can be measured in quantum 
òptics. This required a generalization of the weak value 
theory of Aharonov, Albert and Vaidman. However, in 
a special case of considerable interest, the Stony-Brook 
wave-particle cavity QED experiment, I show that the 
correlation function reduces to a form as simple, and al- 
most the same, as that originally derived by AAV. More- 
over, this experiment measured not a weak value at a 
particular time, but the change in a weak value over time. 
That is, it showed the dynamics of a weak value. 

Weak vàlues are of considerable intrinsic theoretical 
and philosophical interest. Moreover, the fact that they 
obey certain well-defined rules [[|||] hnply that they are 
also heuristic tools for predicting the results of experi- 
ments. Of course the results could be obtained by other, 
more laborious methods. But the power of having several 
fonts of intuition should not be undervalued. The analy- 
sis of Sec. IV is a case in point. This work should help to 
identify other instances in which the theory of weak vàl- 
ues can be used in quantum òptics and related areas. It 
may even suggest new experiments that provide further 
illustrations and applications for this thought-provoking 
area of quantum physics. 
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